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ABSTRACT 

We present results from Chandra/ ACIS-S, Spitzer, XMM-Newton, HST, and VLA 
observations of the radio hot spots, extended environment, and nucleus of the nearby 
(;z=0.0597) FR II radio galaxy 3C 33. This is a relatively low -power FR II radio 
galaxy, and so we expect, a priori, to detect a significant X-ray synchrotron compo- 
nent to the emission from the hot spots. We detect X-ray emission coincident with 
the two knots of polarized optical emission from the southern hot spot (SHS), as well 
as along the northwest arm of this hot spot. We also detect X-ray emission from two 
compact regions of the northern hot spot (NHS), as well as diffuse emission behind 
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the radio peak. The X-ray flux density of the region at the tip of the southern hot spot, 
the most compact radio feature of the southern lobe, is consistent with the synchrotron 
self-Compton (SSC) process. The X-ray flux densities of the other three regions of 
the SHS and the two compact regions of the NHS are an order of magnitude or more 
above the predictions from either the SSC and inverse-Compton scattering of the CMB 
(IC/CMB) mechanisms, thus strongly disfavoring these scenarios unless they are far 
from equipartition (B ~4-14 times smaller than the equipartition values). The X-ray 
flux from the diffuse region behind the NHS is consistent with the IC/CMB prediction 
assuming a small departure from equipartition. However, the radio to X-ray flux den- 
sity distributions of the three regions of the SHS cannot be fitted with a fully convex, 
monotonically decreasing spectral model as a function of frequency. We conclude that 
the X-ray emission is synchrotron emission from multiple populations of ultrarelativis- 
tic electrons unless these regions are far from equipartition. There must therefore be 
complex, unresolved substructure within each knot, similar to that which is found in 
Chandra observations of nearby FR I jets such as Centaurus A and 3C 66B to explain 
the flux density distributions. The nuclear spectrum consists of a heavily absorbed 
(N H ~4x 10 23 cm~ 2 ) power law with a bright Fe K a line (EW ~200 eV) and an Fe 
K edge, thermal emission from galactic gas, and an additional component that could 
either be reflection from the torus or a second power law representing emission from 
an unresolved pc scale jet. The detection of X-ray emission from both hot spots com- 
bined with the large absorbing column toward the primary power-law component of 
the nucleus conclusively demonstrate that the jets must lie relatively close to the plane 
of the sky and that relativistic beaming cannot be important. We also detect diffuse 
X-ray emission, some of which is likely to originate in hot galaxy or group-scale gas, 
and some of which may be attributed to IC/CMB off the relativistic electrons of the 
radio lobes. 

Subject headings: galaxies: individual (3C33) - X-rays: galaxies: clusters - galaxies: 
ISM - hydrodynamics - galaxies: jets 



1. Introduction 



The origin of X-ray emission from the radio hot spots of FR II radio galaxies is still a subject 
of considerable debate. For many sources, inverse-Compton scattering of radio photons by the 
relativistic electrons of the hot spot (the synchrotron self-Compton mechanism (SSC)) with the 
magnetic field at or near the equipartition value provides an adequate description o f the broad 



band (radio to X-ray) spectra (|Harris. Carilli. & Perley II 19941 : lHardcastle et al. 2002). There are 
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however, many sources for which the X-ray emission is orders of magnitude larger than the SSC 
prediction and for which this paradigm is clearly incomplete. It has recently been suggested, based 
on analysis of archival Chandra observations of the compact hot spots of 3CRR radio galaxies, 
that in some cases synchrotron emission c an be the dominant co mponent and can be orders of 
magnitude larger than the SSC component (IHardcastle et al. Il2004|) . It was found t hat lower radio 



power hot spots tended t o have stronger optical and X-ray synchrotron components (|Brunetti et al. 



20031 : IHardcastle et al. Il2004|) . It was speculated that this represents an intermediate case between 
X-ray synchrotron emission as seen in the jets of FR I radio galaxies and the SSC mechanism 
of the powerful FR II hot spots. Some hot spots do indeed appear to show a SEP consistent 



with synchrotron emission from a simple continu ous -injection mod el (|Heavens & Meisenheimer 



9871) between the radio and X-ray (e.g. 3C 403: iKraft et al.ll2005|) . However. iTavecchio et al 



(|2005|) showed that good optical/UV constraints on the fluxes of the hot spots in some FR II quasars 
preclude fitting a single-component synchrotron model to the radio through X-ray data, while the 
normalization of the X-ray emission is inconsistent (orders of magnitude too large) with an SSC 
model. Instead, they suggested th at models involving beaming (such as the deceleration model of 
Georganopoulos & Kazanasll2003|) may be required to explain the observations. The synchrotron 
hypothesis can be directly evaluated with high resolution X-ray images of the hot spots of nearby 
FR lis: high spatial resolution is necessary if we are to understand the relationship between the X- 
ray emission and the various radio and optical structures related to jet termination in these sources. 

In the paper, we present results from new Chandra, Spitzer, and XMM-Newton observations, 
combined with archival HST, GALEX, and VLA observ ations, of the hot spo ts, the nucleus, and 
the extended environment of the nearby (,2=0.0597, see (|Popescu et al. II 1992 )) FR II radio galaxy 
3C 33. 3C 33 is a double radio source that extends ~4' on the sky with a weak radio core 
( Hargrave & McEllin Il975). The radio emission at 1 .4 GHz is dominated by the southern hot 
spot (SHS) (|Rudnick 1988 : Rudnick & Anderson 1990 ). Optical emission coincident with the ra- 
dio peak of the SHS was first reported by Isimkinl (|l979 ). and it was subsequently determined that 
this emission is highly polarized (~40%), thus conclusi vely confirming it as synchrotron emission 
(|Meisenheimer & Roserlll986l : ICrane & Stiavelli Ill992[) . Imaging and spectroscopi c studies of the 
hot spot have shown a complex o ptical morphology with a steep power-law slope (|Simkin 1 ll986l: 



Crane. Stockton. & Saslow 



1987|) . The hot spots were not detected in a ROSAT/HRI observation, 



but the upper limit to the X-ray flux demonstr ated that there must be a break in the sy nchrotron 
spectrum between the optical and X-ray bands (IHardcastle. Birkinshaw. & Worrallll 19981) . The pri- 
mary scientific goal of the new observations is to resolve spatially the emission regions of the hot 
spots from the radio to the X-ray band in order to elucidate the emission processes that are im- 
portant at the terminal hot spots of FR lis. 3C 33 is an ideal candidate for this study as it is the 
nearest nar row-line FR II in the 3CRR c atalog with optical synchrotron emission detected from 
a hot spot (|Meisenheimer & RoserHl986|) . The proximity of 3C 33 permits us to resolve spatially 
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the emission features at the highest possible linear resolution. In addition, since this source has a 
narrow-line nucleus, relativistic beaming is not likely to be important in modifying the spectral or 
morphological properties of the hot spots. 

This paper is organized as follows. Section 2 contains a summary of the observations, the 
results of the data analysis are presented in section 3, and we discuss the implications of our 
results in section 4. Section 5 contains a brief summary an d conclusions. We assume WMAP 
Year-1 cosmology throughout this paper (|Spergel et al. 112003b . The observed redshift (2=0.0597) 
of the host galaxy of 3C 33 corresponds to a luminosity distance of 263.9 Mpc, and one arcsecond 
is 1.14 kpc. All uncertainties are at 90% confidence for one parameter of interest unless otherwise 
stated, and all coordinates are J2000. Absorption by gas in our Galaxy (iV#=4.0x 10 20 cm -2 ) 



(|Dickey & Lockman Ill990() is included in all spectral fits and count rate to flux density conversions. 



2. Observations 

The radio galaxy 3C 33 was observed with Chandra/ ACIS-S in two observations (OBSID 
6190 and 7200) on November 8, 2005 and November 12, 2005 (PI: Stephen Murray). The total 
observing time was ~40.3 ks. The observation was pointed so that the SHS was near the best 
focus, but all regions of the radio galaxy, including the northern hot spot (NHS) and nucleus, 
were contained on the S3 chip. We made light curves for the S3 CCD in the 5.0 to 10.0 keV 
band, excluding the nucleus, to search for background flares, and removed intervals where the 
background rate was more than 3a above the mean, leaving 39831.6 s of good data. Bad pixels, 
hot columns, and events along node boundaries also were removed, and standard ASCA grade 
filtering (0,2,3,4,6) was applied to the events file. Emission from 3C 33 fills only a small fraction 
of the field of view of the S3 chip, so local background was used in all spectral analysis. The 
X-ray core is aligned to ~0.1" with the radio core, so no adjustment of the X-ray coordinates was 
necessary. All X-ray images presented in this paper were generated from the Chandra data. 

XMM-Newton observed 3C 33 for a total of 17 ks in two separate observations on 2004 Jan- 
uary 4 and 2004 January 21 (PI: Martin Hardcastle). The data were processed with the Scientific 
Analysis Software (SAS) using the standard pipeline tasks emchain and epchain and filtered for 
patterns <= 4 (MOS) or <= 12 (pn) and using the bit-mask flags 0x766a0600 for MOS and 
OxfaOOOc for pn, which are equivalent to the standard flagset #XMMEA_EM/EP but include out of 
field-of-view events and exclude bad columns and rows. In the first of the two observations, the 
background count rate was low (~0.2 counts s _1 for MOS and ~0.4 counts s _1 for pn) so that no 
temporal filtering was required; however, both the MOS and pn data from the second observation 
suffered from a large flare at the start of the observation, which rendered the pn data unusable and 
required the MOS data to be filtered using a count rate threshold of 0.5 counts s _1 . The resulting 
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clean, merged data had a duration of 9700 s (MOS) and 6314 s (pn). The XMM-Newton data were 
used primarily in our study of the nucleus and the extended, low surface brightness diffuse emis- 
sion from the group. The SHS was detected in these observations, the NHS was not. The larger 
XMM-Newton PSF (relative to Chandra) provides no additional structural constraints, so we use 
only the Chandra data in our study of the hot spots. 

We present data from three Spitzer observations of 3C 33, two data sets taken from the public 
archive (one IRAC and one MIPS observation, both pointed towards the SHS - program ID 3327), 
and an additional IRAC observation (program ID 3418) centered on the host galaxy taken by the 
3CRR low-z consortium (Birkinshaw et al. 2007, in preparation). The IRAC observations were 
made on July 23, 2004 (3327), and January 6, 2005 (3418), with observations times of 96 and 
360 seconds, respectively. The MIPS (24 fim only) observation was made on December 23, 2004 
(3327) with an observation time of 200 seconds. The images used were produced by the data 
analysis pipeline version 14.1.0. The SHS is detected in all four IRAC bands and the MIPS image. 
The NHS is only contained within the Spitzer field of view in the 2005 IRAC observation and the 
MIPS image, and is detected in all four IRAC bands and the MIPS 24 fxm band, although there is 
some confusion from adjacent stars. 

We also use archival HST (WFPC and WFPC2) and VLA data in this paper. The SHS 
was observed for 2400 s with the //ST/WFPC2 instrument in 1995 as part of the HST survey 
of radio-galaxy hot spots (PI: P. Crane) using the F702W filter (pivot wavelength of 6919 A), 
and for 1 800 s with the HSTfWF FC instrument using the F606W filter (pivot wavelength of 



5888 A) (ICrane & Stiavelli Ill992h . There have been no HST observations of the NHS. We ob- 
tained the reprocessed data of the SHS from the HST archive and used the IRAF synphot package 
to apply photometric calibrations. T he fluxes were reddening-corrected using the dust maps of 



Schlegel. Finkbeiner. & Davisl dl998h . The correction to the flux densities is ~25% in the visual. 



A 1.5-GHz r adio map with a r esolut ion of 4.0 arcsec was obtained from the 3CRR Atlas^: this 
the image of Leahy & Perleyl (|l99ll) . At higher frequencies, we used 4.9 and 15-GHz data from 



the VLA public archive. Details of the VLA data used and the maps made from them are given in 
Table [TJ The data were calibrated and reduced in the st andard manner usin g AIPS. Results from 



analys es of these data have been previously published in lRudnickl (|1988|) and lRudnick & Anderson 



(11990) . 



Finally, we also examined the XMM-Newton Optical Monitor (UVM2 filter) data and archival 
GALEX images of 3C 33. The hot spots were not detected in any of the UV images. The 
observation time of the GALEX observation was sufficiently short (~90 s) that the upper lim- 
its are not interesting (they lie well above an extrapolation of the optical detections or limits). 



1 http : // w w w. jb . man . ac . uk/atlas/ 
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The XMM-Newton OM upper limits lie above the GALEX limits because of a combination of 
the narrower bandpass and lower effective area. The 3a upper limits on the flux density of 
SI in the GALEX NUV (pivot wavelength 2297 A) and FUV (pivot wavelength 1524 A) bands 
are 2.0 and 6.1 aJv, respect i vely. These values were dereddened using the extinction curve of 



Cardelli. Clavton. & Mathis (1989). 



3. Analysis 
3.1. Hot Spots 

An adaptively- smoothed Chandra/ ACIS-S image of 3C 33 with 1.4 GHz radio contours over- 
laid is shown in Figure [TJ The positions of the nucleus and the two hot spots (NHS and SHS) are 
labeled. Unsmoothed Chandra/ ACIS-S images of the SHS and NHS in the 0.5-5.0 keV band with 
4.9 GHz radio contours overlaid are shown in Figures [2] and [3l respectively. We have divided the 
SHS and the NHS into four and three regions, respectively, for analysis. The regions S1-S4 and 
Nl-2 are shown in the two figures and summarized in Tabled Two additional regions, labeled 
DN1 and DS1, coincident with diffuse X-ray emission between the nucleus and the radio hot spots 
are also contained in Table[2] The 1 keV flux density is given for each region assuming a power law 
spectrum with photon index 2 and Galactic absorption, consistent with the spectral fits described 
below. 

In both the NHS and SHS, spatially-resolved X-ray emission is detected coincident with the 
compact radio peaks (regions SI and Nl). In the SHS, X-ray emission is also detected closer 
to the nucleus than the hot spot to the northeast (region S2), and along the entire northwest half 
of the lobe (regions S3 and S4). An #ST/WFPC2 image of the SHS with 4.9 GHz radio contours 
overlaid is shown in FigureHl An optical peak coincident with the radio peak (S 1) is observed, as is 
a weaker feature ~2" to the northeast (S2). Both of these features have X-ray counterparts. There 
is no apparent optical emission from S3 or S4. In the NHS, a second region of X-ray emission is 
detected in the NHS to the northwest of the radio peak (labeled region N2). There is also a tail of 
diffuse X-ray emission behind (i.e. toward the nucleus) Nl and N2. This can be seen in Figure Q] 
and has been labeled N3 in Table [2] 

Before proceeding, we must clarify an important point of terminology. Historically the bright 
radio features of 3C 33 (in their entirety) that we have termed the NHS and SHS in this paper 
have been generally been referred to as 'hot spots' because they are relatively compact. This 
nomenclature is not entirely correct, however, as they contain both compact features and resolved, 
diffuse emission behind the compact features. Neither are they 'lobes', however, such as are seen 
in the canonical FR lis Cyg A and 3C 98. For historical continuity, we will continue to use the 
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acronyms NHS and SHS or the term 'hot sp ot' to refer to the en tire radio bright features, but we 



will use the definitions and nomenclature of [Leahy et al. I (| 19971) when discussing our individual 



regions. By these criteria, SI and Nl are primary hot spots. They are compact and are reasonably 
associated with the termination of the jet. Regions S2 and N2 are secondary hot spots. They are 
clearly hot spots in a morphological sense, although they are not the termination points of the jet. 
Finally, regions S3, S4, and N3 are considered to be part of the 'hot spot complex', although they 
are not hotspots. We will also use the term 'hot spot' to refer to compact features of FR lis and 
radio quasars that are believed to be the termination of the jet. The term 'hot spot' will not be used 
in connection with individual regions that are clearly not compact (e.g. S3, S4, and N3). 

There are too few counts to fit the X-ray spectra of each region of the NHS and the SHS 
individually (typically 10-15 source counts in each region), so we combined all the counts in Sl- 
S4 and N1-N3 to create two spectra. We fit each spectrum with an absorbed power law model, 
with the absorption fixed at the Galactic value to reduce the number of free parameters. The best 
fit photon indices are 1.8±0.6 and 2.2±0.8 for the SHS and NHS, respectively. The quality of 
these data are not sufficient to discriminate between emission models (i.e. synchrotron, thermal, 
or inverse Compton scattering) on the basis of spectra, although we reject the thermal model for 
reasons detailed below. The unabsorbed X-ray luminosities of the SHS and NHS are ~1.6 and 
~0.5x 10 41 ergs s -1 , respectively, in the 0.5 to 5.0 keV bandpass. 

The morphology of the SHS in the infra-red in two Spitzer/IRAC bands, 4.5 and 8.0 /im, is 
shown in Figures [5] At the shortest IRAC wavelengths (Figure [5] left), region SI dominates the 
IR emission of the SHS, although S2 is clearly detected. There is little or no emission from S3 or 
S4. At longer wavelengths (Figure \5\ right), the fluxes from SI and S2 are roughly equal. Thus the 
flux density spectral index (and therefore the electron energy distribution index) is considerably 
steeper for S2 than S 1 . The morphology of the NHS in the IR at 4.5 /im, shown in Figured is very 
different than that of the SHS. There is a bright IR peak coincident with region Nl, and a ~8" long 
tail of emission to the west that follows the radio contours through region N2. The limited spatial 
resolution of Spitzer (~1.2") prevents us from determining whether this emission is continuous or 
composed of small (kpc-scale) knots. If the radio peak (Nl) represents the termination point of 
the jet, the detection of this filament of IR emission perpendicular to the hot spot/jet axis suggests 
that there is particle acceleration (and therefore large velocities) perpendicular to this axis as well. 
This is quali tatively consistent with the 'p rimary + splash model' for the flow pattern of the NHS 
proposed by Rudnick & Anderson] (|l990|) . 



The angular resolution of Spitzer is such that the SI and S2 components of the SHS, and the 
Nl and N2 components of the NHS, are separated in IRAC channels 1 and 2 (3.6 and 4.5 /im), but 
blended to a greater or lesser extent at the other wavelengths. The fields are further complicated by 
the presence of nearby foreground stars. To extract the flux densities of the separate components 
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we therefore fitted point sources at the approximate locations of the components (leaving some 
freedom in the coordinates to allow for a misalignment of the radio and IR coordinate systems), 
and at one or two contaminating stars in the hot-spot fields. Cross-checks of the derived flux 
densities were made using aperture photometry, where possible. Errors on the flux densities of 
the components were derived from the fits taking account of the residuals in the areas fitted and 
the number of parameters being fitted, and include the systematic errors on the flux density scale 
expected in the different bands. 

At 24/mi this degree of fitting was not possible, because the hot spot components are severely 
blended by the diffraction-limited point spread of Spitzer. At this band we therefore fitted the best 
flux densities using defined locations for the hot spot components, and compared the results with 
the fraction of the total flux densities for the hot spot regions that would be expected given their 
3.6 to 8 /im spectra. Generally the results agree, though the errors from the spectral extrapolation 
constrained by the total flux densities of the hot spots are smaller, and so the latter results were 
adopted. 

The results for the IR flux densities, as well as radio flux densities, of all seven regions (SI 
through S4 and Nl through N3) are shown in Table [3] For completeness, the optical, UV, and 
X-ray flux densities (or upper limits) for SI are shown in Table [4] The optical points are taken 



from either archival HS T images or previously published values (|Meisenheimer & Rosen 1 19861 ; 



Crane & StiaveUilll992|) . Regions S3, S4, and Nl through N3 were not detected in any optical or 



UV band. 

We have constructed radio through X-ray flux density distributions for six regions. These are 
shown in Figures [7J [8l [9] for the Sl/2, S3/4, and Nl/2 pairs, respectively. All error bars are la 
statistical uncertainties. If the error bars are not visible, they are smaller than the data points. The 
upper limits on the IRAC fluxes shown in Figure[8]are 3a limits. The solid lines are representative 
single injection synchrotron models, and the dotted and dashed lines are the estimates of the mini- 
mum energy IC/CMB and SSC models, respectively. The details of these spectral models, as well 
as implications regarding X-ray emission mechanisms, are presented in Section 4. 



3.2. Nucleus 

X-ray images of the nuclear region of 3C 33 in the 0.3-0.7 (soft) and 1.5-2.0 (hard) keV bands 
are shown in Figure QTJ In the hard band, the nucleus is unresolved and consistent with a point 
source. In the soft band the structure is clearly extended, with the emission elongated along the 
northeast/southwest axis. This extension is not the frame transfer streak (FTS) as it lies roughly 
along the ENE/WSW axis through the nucleus. The position angle of the observed extension is 
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rotated relative to this axis by approximately 30°. In addition, the nucleus contains only ~2000 
counts, so ~30 counts are expected over the entire ~8' length of the FTS. HST observations of 
the host galaxy show a similar elongation in the stellar distribution suggesting that the extended 
X-ray emission is from thermal gas. There is an additional spatially extended X-ray component 
~5" (5.7 kpc) to the northeast. It is unclear whether this is thermal emission from an extended 
gaseous component or X-ray emission from a nuclear jet. There is no evidence of a kpc scale jet 
in any of our radio maps. 

The spectrum of the nucleus was extracted from a source-centered circle of radius 10 pixels 
(4.92"), with background sampled from a surrounding annulus of inner radius 10 pixels and outer 
radius 20 pixels (9.84"). The Chandra count rate from the nucleus is ~0.05 cts s~\ indicating 
little pileup. Initially we fitted the spectrum of the nucleus with a three-component model: an 
absorbed power law, a Gaussian line representing emission from neutral Fe K a , and a thermal 
(the XSPEC Astrophysical Plasma Emission Code model, or simply APEC) component. The fit 
was poor (xl=l9.9 for 26 degrees of freedom) with large residuals between 1 and 3 keV, so we 
investigated two alternative models. 

Model 1 consisted of a four component model, the three components listed in the previous 
paragraph plus a reflection component. The reflection model has been used to fit the X-ray spec- 
tra of Seyfert galaxies, but has rarely been found to provide an adequate description of the X-ray 
spectra of radio-galaxy nuclei. The X-ray spectrum with the best-fit four-component reflection and 
two power-law component models and residuals are shown in the left and right plots of Figure [T4l 
respectively. Model 2 also consisted of the three components listed in the previous paragraph 
plus a second, less absorbed power law representing emission from a parsec scale jet. Such a 
model has bee n shown to be a good description of emission from the nuclei of FR II radio galaxies 



(lEvans et al. II2006Q . A summary of best-fit parameters and uncertainties for both models is con- 
tained in Table[5] The APEC component in the model that includes a reflection component required 
abundances of Mg and Si that were larger than the Solar value and a sub-Solar value for Fe. This is 
probably not physical and may be indicative of systematics in the model such as multi-temperature 
gas in the ISM or more complex reflection geometry than a uniformly illuminated slab. We note 
that the large column density toward the primary power law in both spectral fits is consistent with 
the expectations from unified models for a narrow-li ne radio galaxy, as is generally seen in X-ray 



observations of large samples of such so urces (e.g. ISambruna et al. Ill999i : LBelsole et al. 1120061 : 



Evans et al. 1 12006: Ha rdcastle et al. l bood) . This strongly supports the idea that the jet axis is close 



to the plane of the sky and that relativistic beaming is not important in modifying the appearance 
of the hot spots. Finally, we mention that we fitted the XMM-Newton spectra (both MOS cameras 
and PN) with the same models and find that the fits are consistent with the Chandra results. 



- 10- 



3.3. Large Scale Diffuse Emission and Environment 

Diffuse, low surface brightness emission unrelated to either the hot spots, the central galaxy, 
or the active nucleus is visible in Figure [TJ This emission is more clearly shown in Figure [T2l an 
adaptively smoothed X-ray image in the 0.5-2.0 keV band with all point sources other than the 
NHS and SHS removed. Two regions of diffuse X-ray emission have been labeled DN1 and DS1. 
This emission is generally aligned along the NE/SW axis of the radio structure, and extends more 
than 100 kpc. There is a rough spatial correspondence between this X-ray emission and the diffuse 
radio lobe emission between the hot spots and the nucleus. The X-ray fluxes of regions DN1 and 
DS1 are 4.5 and 6.2 xlCT 15 ergs cm~ 2 s _1 assuming a power-law spectrum with photon index 1.5 
and Galactic absorption: this corresponds to 1-keV flux densities of 1.3 and 1.8 nJy. The combined 
spectrum of the two regions is well fitted with a power-law model of photon index 1.5±0.7 (90% 
confidence) and Galactic absorption. We also fitted the spectrum with an absorbed APEC model 
with the absorption fixed at the Galactic value and the elemental abundance fixed to 0.5 times Solar 
to reduce the number of free parameters. The best fit temperature is l.Sljjjjj keV (90% confidence 
for one parameter of interest). We believe that this emission is non-thermal, but a purely thermal 
model is not conclusively ruled out by the data (see Section 4.2 for complete discussion). 

To investigate the properties of the undisturbed environment of 3C 33 on smaller (i.e. galactic) 
scales, we extracted a radial surface brightness profile from the Chandra data in the 0.5 - 2 keV 
energy range, excluding the radio lobe and hotspot regions. This profile is shown in Figure[[3l The 
total number of counts (after background subtraction) in the surface brightness profile is 140^28 
cts. We modeled the Chandra point-spread function using ChaRT and Marx. A point-source 
model was not a good fit to the profile, with a clear detection of extended emission between 10" 
and ~100". We therefore fitted a point-source plus [3 model to the profile, although the parameters 
of the (3 model were poorly constrained. Acceptable fits were obtained for (3 values in the range 
of 0.67 to 0.9, with r c ~2" (~2.3 kpc). We extracted a spectrum in an annulus between 10 and 
100 arcsec, excluding the radio lobes and hotspots, and found that the 0.5 - 2 keV spectrum is 
consistent with a thermal model of /csT=1.8^o g keV (the spectrum above 2 keV is dominated by 
scattered AGN emission). This temperature is unphysically large for a relatively isolated elliptical 
galaxy and perhaps suggests that even the gas in the central 10 kpc of the host galaxy has been 
disturbed by the passage of the powerful radio jets. We used the XMM-Newton data to investigate 
whether any emission on larger scales was present; the data are consistent with the presence of 
the Chandra-detected component, but the poor data quality meant that no strong upper limits on 
larger-scale emission could be obtained. 

We integrated the best-fitting (3 model {(3 = 0.9, r c = 1.3") to determine the total net counts 
from the environment, which corresponds to a count rate of ~7x 10~ 3 cts s _1 in the 0.5-2.0 keV 
band. Using the best-fitting temperature, this corresponds to a luminosity of 1.5 x 10 41 ergs 
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s 1 in the 0.5-2.0 keV band, or a bolometric lu minosity of ~3xl0 41 ergs s" 1 . The lum inosity- 
temperature relation for radio-quiet groups from lCroston. Hardcastle. Birkinshaw I (120051) predicts 
a luminosity of ~4xl0 43 ergs s _1 for groups with a temperature of ~1.8 keV, so that the 3C33 
group appears to be roughly two orders of magnitude less luminous than would be expected for 
its temperature, unless both Chandra and XMM-Newton are failing to detect a very large, smooth 
extended environment. A temperature of ~0.5 keV is predicted for our measured luminosity. This 
suggests either that 3C 33 (or previous generations of outbursts from its AGN) has had a dramatic 
effect on the surrounding group, or that powerful radio galaxies can form in extremely poor envi- 
ronments. 



4. Interpretation 

4.1. X-ray Emission from Hot Spots 

We have fitted three models to the radio through X-ray flux density distributions of the seven 
regions shown in Figures [7J through [lOl a broken power-law synchrotron model with a high en- 
ergy cut-off, a synchrotron self-Compton model (SSC) in which the X-ray emission is the result 
of inverse-Compton scattering of the synchrotron radio photons off the relativistic electrons, and a 
model of inverse Compton scattering of CMB photons (IC/CMB) off the radio synchrotron emit- 
ting relativistic electrons. Other sources of seed photons in the inverse-Compton scenario including 
optical photons from the stellar population of the host galaxy or beamed optical/IR photons from 
a hidden nuclear jet are implausible beca use of the distance of the hot spots from the nucleus (see 



Hardcastle. Birkinshaw. & Worra ll (1998) for details). For example, the energy density of IR pho- 
tons emitted from the nucleus (based on the nuclear flux measured in the Spitzer images) is roughly 
two orders of magnitude below the energy density of the CMB at the distance of the hot spots. 

The relevant model parameters for the synchrotron, IC/CMB, and SSC curves shown in Fig- 
ures [7J through [TO] are summarized in Table|6j We assume a power-law electron energy distribution 
for the synchrotron fits with 7 mi „=100 and 7 ma a;=10 8 . We also allow for a break in the electron 
energy distribution where the index of the power law steepens. We choose to plot synchrotron 
models for S3 and S4 with relatively large breaks for consistency with SI and S2. For the latter 
two regions, the large change in power-law index is required by the optical and IR data. The flux 
density distributions of S3 and S4 could be modeled with a smaller change in index which would 
be roughly consistent with the optical and IR limits, although the lack of detections in these bands 
prevent a definitive conclusion. We emphasize that in regions S3 and S4, we must be averaging 
over considerable substructure in the synchrotron scenario because of the short lifetime of the par- 
ticles. The broken power-law model thus should be viewed as a parameterization of the data and 
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does not imply a single-zone model. For the IC/CMB and SSC fits, we assume equipartition with 
no contribution from relativistic protons. For all regions other than SI, we assume a cylindrical 
geometry in the plane of the sky with the diameter of the cylinder equal to the short axis of the 
region, and the height equal to the long axis. The IC/CMB fits roll off at low frequency because of 
the choice of 7 m j n . 

The SSC calculation is sensitive to the assumed geometry of the spot, and high-frequency ra- 
dio observations of S 1 clearly indicate the presence of compact structures smaller than the box cho- 
sen for X-ray analysis. We consider two models of the geometry of S 1 for SSC calculations. Firstly, 
we assumed that all the radio emission in the SI region comes from a uniform-density sphere of 
radius Of! 5. This value was chosen as a reasonably compromise be tween the the 15 GHz radio mor- 



phol ogy and the much smaller value of Of! 1 7 used in the analysis qflHardcastle. Birkinshaw. & Worrall 



(I1998|) . who based it on the modeling of iMeisenheimer et al. I (|1989|) . The results of this simple 



model are plotted in Figure [71 As the hotspot region is clearly resolved in both radio and X-ray, we 



next used the code of lHardcastle et al. I (|2002|) to consider a more complex two-zone model which 
better matches the structure of the hotspot in the highest-resolution 15-GHz observations. This 
model consisted of a section of a thin spherical shell with uniform electron density embedded in a 
region bounded by a paraboloid with a linear electron density gradient: the magnetic field strengths 
in the two components were separately normalized using their observed 15-GHz flux densities, and 
other electron spectral parameters were as described above. The predicted X-ray flux density of 
this model (which is dominated by the more extended filled paraboloid region, as a thin shell is a 
poor source of SSC emission) is roughly a factor 2 below the estimate derived from the spherical 
geometry. We conclude that the largest uncertainty in the SSC flux from SI is due to the assumed 
geometry, but that the observed X-ray flux is consistent with the SSC prediction with only a small 
departure from equipartition. 

Although the observed X-ray flux from SI is consistent with the SSC prediction, the X-ray 
flux densities of regions S2, S3, S4, Nl, and N2 are one to two orders of magnitude larger than 
predicted in either the SSC and IC/CMB scenarios if the features are in equipartition. (see Figures[7J 
through [TOl). The larger region N3 is reasonably well described by the IC/CMB model with a small 
departure from equipartition. The radio through optical flux density distributions of the regions S2 
through S4 indicate a significant steepening of the spectral index in the optical region, so that the 
observed X-ray emission cannot be described by an extrapolation of a simple single-component 
synchrotron model for these regions of the SHS; if this extrapolation is plotted (e.g. in Fig. [7J it 
would clearly violate the optical constraints on the flux density distributions. In fact, it has been 
argued that this radio to optical spec tral steepening strengthens the case for a high-energy cut-off 



in the electron energy distribution (|Meisenheimer et al. Ill989i : lHardcastle. Birkinshaw. & Worrall 



19981). Thus, while the X-ray emission from SI is consistent with the SSC mechanism, neither of 



the two canonical emission mechanisms provide a satisfactory explanation for the X-ray emission 
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observed in the other regions of the SHS of 3C 33. The X-ray emission from the two compact 
regions of the NHS (Nl and N2) is orders of magnitude above the SSC prediction, but is well 
described by a simple single-component synchrotron model (shown on Figure [9]). 

There are several important implicit assumptions in the various IC scenarios that may not be 
appropriate for the regions of 3C 33 where SSC and IC/CMB are nominally inadequate. Such fail- 
ures could account for the discrepancies between the observed and predicted X-ray flux densities. 
In particular, it is possible that these regions are not in equipartition, that protons make a signif- 
icant contribution to the energy density, or that the filling factor of the relativistic plasma is less 
than unity. Only a comparatively small departure from equipartition, a factor of ~4 in magnetic 
field (in the sense that B « £> cq /4), is required for regions S2 through S4 to be fitted with an 
IC model. However, the departure is considerably larger for the two compact regions of the NHS 
(B ?s S eq /14). The equipartition expectations are only reduced if there is an energetically signif- 
icant proton population. In principle arbitrarily large SSC fluxes may be produced by low filling 
factors, but in practice these would have to be implausibly low, with volume filling factor values 
~ 10~ 9 , to explain the two orders of magnitude discrepancy between the predicti on and the obser- 



vation s in components such as S4. These issues are explained in more detail in lHardcastle et al. 
d2004l) : lKataoka & Stawarz1d2005l) . 



Relativistic beaming is often invoked to explain the X-ray emis sion of jets of radio quasars 



and t he hot spots of FR II radio galaxies in the various IC scenarios (|Georganopoulos & Kazanas 



2003|) . This is untenable in regions S2-S4 and N1-N2 of 3C 33 for at least five reasons. First, and 
most importantly, X-ray emission has been detected from both the NHS and SHS. It is unlikely 
that both are moving towards us at relativistic velocity near our line of sight. Second, the X- 
ray nucleus is heavily obscured, suggesting that the jets lie close to the plane of the sky. Third, 
the optical nucleus is classified as a narrow emission line radio galaxy, supporting the previous 
point. Fourth, core prominence (the ratio of 5 GHz luminosity of the core to 178 MHz flux of 
the entire source) is often used to estimate th e role of relativistic beaming (lOrr & Brownelll982 : 
Morganti et al. 111997c lHardcastle et al. Il2004|) . The core prominence of the 3C 33 nucleus is low, 
supporting the conclusion that beaming is not important. Fifth, VLBA observations of the core 
of 3C 33 show only a small (factor of 2 at a distance of 5 mas from the core with the brighter 
compon ent aligned with the SH S) asymmetry in the ratio of the flux density of the jet/counterjet at 
5 GHz (|Giovannini et al. Il2005|) . Thus, we conclude that relativistic beaming is not an important 
factor. 

As an alternative, we consider a scenario in which internal motions in the regions at relativis- 
tic velocities could account for the observed X-ray emission. In this model, the radio emission 
originates in two distinct components, one which follows the general flow of the hot spots (i.e. 
non-relativistic in or near the plane of the sky) and a second component that has been deflected 
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close to our line of sight at relativistic velocity. To evaluate the plausibility of this model, we con- 
sider knot S4. If we assume that all of the radio emission from this region originates in the beamed 
component and it occupies the entire volume of S4, a bulk Lorentz factor, T, of ~2 or greater is 
required, and the direction of motion must be close to the line of sight (9 ~ 1/T). Assuming a 
smaller volume would only push T higher. There is some evidence for mildly relativistic flows 
behind FR II hot spots, although nothing of this magnitude. We conclude that it is unlikely that the 
X-ray emission from these regions can be accounted for by relativistic motions of sub-components. 

It is also extremely unlikely that the X-ray emission from these regions is due to hot gas that 
has been swept up, since the host galaxy of 3C 33 resides in a poor environment. Assuming that the 
X-ray emission of the SHS is from a thermal plasma with k B T=l .5 keV and Z=0.5 times Solar, the 
density of the gas would be ~10 -2 cm~ 3 assuming a uniform filling factor. The pressure and mass 
of this gas would be ~6.2x 10~ n dyn cm~ 2 and ~10 9 M & , respectively. The gas pressure of this 
putative shell would e xceed the equi partition pressure of the diffuse regions behind the hot spot 



(SI) by a factor of ~6 (|Rudnicklll988l) . This is not implausibly large as the equipartition pressure 



of lobes and jets of FR I radio galax ies is often less th an the pressure in the ambient medium (e.g. 



the inner radio lobes of Cen A, see lKraft et al. 1120031) . although it would be unusual for an FR II 



If the extended diffuse X-ray emission described in subsection 3.3 is from an extended corona, the 
mass of gas in the NHS and SHS would each be ~5% of the total gas mass, which is implausibly 
large. 

The only other possibility is that the X-ray emission is synchrotron radiation from a pop- 
ulation of ultra-relativistic electrons, although the observed flux density distributions cannot be 
simply described with single or continuous injection models. There are, in fact, several arguments 
that support the multi-zone synchrotron hypothesis for regions S2, S3, S4, Nl, and N2. First, 
the radiative lifetime of the ultra-relativistic particles is a few tens of years in the equipartition 
magnetic fields of the regions. All of the knots are resolved in the X-rays, so the particles must be 
re-energized along the length of the regions. Thus, in each knot there must be many sites of particle 
acceleration. There is therefore no reason to expect that a single injection model will accurately 
describe the broad band flux density distribution from such a scenario. Second, it is clear from 
Chandra observations of nearby FR I jets, whose radio through X-ray emission is believed to be 
synchrotron, that the morphological relation ships between the radio, IR, optical, and X-ray syn 



chrotron emitting electrons are complicated ([Kraft et al]|200l] ; lHardcastle. Birkinshaw. & Worrall 



2001). 



It is clear, however, that the simple model described in Har dcastle et al. I (|2004), in which a 
single continuous-injection spectrum describes the radio through X-ray spectrum of these regions, 
does not apply to 3C 33. The resolved nature of the X-ray emission, coming from throughout the 
bright 10-kpc SHS, in fact makes it clear that a model in which the X-rays are generated by elec- 
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trons accelerated at a single location is simply not appropriate in this object. At the same time, both 
the resolved nature of the X-ray emission and the fact t hat 3C 33 lies in the pla ne of the sky rule out 
explanations involving beaming, such as those used by lTavecchio et al. I (|2005|) for their quasar hot 
spots. Given the known complex distributed particle acceleration properties of low-power radio 
jets, we feel it is reasonable to suggest that the SHS of 3C 33 contains many discrete regions of 
particle acceleration and that incompatibility with a one-zone continuous injection model is not in- 
consistent with the X-rays having a synchrotron origin. I n fact , a multiple component synchrotron 
model for the SHS was required by lMeisenheimer et al. I(|1997|) to explain the radio to optical spec- 
tral index and the optical morphology. Deep, high spatial resolution observations in the optical and 
the near-UV may offer the best prospect of testing this picture. 

Therefore, both the synchrotron and SSC processes are important in the SHS of 3C 33. The 
equipartition magnetic field strength is ~195 /iG for region SI (and the field strength would be 
larger for a smaller assumed volume), and between 25-80 fiG for regions S2 through S4. Thus, 
the rate of synchrotron energy loss in region SI is 5-64 times greater than in the rest of the re- 
gions. There is no X-r ay synchrotron emission from SI because the magnetic field is too large 
(|Hardcastle et al. 1120040 . Extended regions of X-ray emission, where the SSC and IC mechanisms 
are out of the question, are strong evidence for a distributed particle acceleration process like that 
seen in low-power jets. This suggests that a simple choice between an SSC or an IC origin for the 
emission mechanism may not be appropriate for other, more distant, hot spots for which we do not 
have the high linear resolution that we have in 3C 33. 



4.2. Diffuse Emission Co-Incident with Radio Lobes 

We consider two possible models for the origin of the diffuse X-ray emission coincident with 
the radio lobes: inverse-Compton scattering of CMB photons off the relativistic electrons in the 
diffuse radio features behind the hot spots and thermal emission from hot gas. The quality of the 
data is not good enough for us to differentiate with certainty between these possibilities. 



4.2.1. IC/CMB 



X-ray emission from the extended lobes of radio galaxies is generally attributed to inverse- 
Compton scatteri ng of cosmic microwave background (CM B) photons by the low-energy elec- 
trons in the lobe (fFeigelson et al. 111995c LKaneda et al. I Il995|) . Additionally, it has been shown 
for large samples of FR II radio sources that the X-ray fluxes measured from the lobes gener- 
ally lie within a factor of a few of the prediction from the inverse-Compton process assuming 
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equipartition magnetic fields ( Kataoka & Stawarzll2005 : Croston et al. I l2005h . We calculate the ex- 
pected inverse-Compton fluxes from the lo bes on the assumption of equiparti tion using the code of 



Hardcastle. Birkinshaw. & Worralll (|1998|) and with the same assumptions as ICroston et al.1 (I2005|) 
(lobes are treated as uniform cylinders, 7 min = 10, 5 = 2, and a break in the electron spectrum is 
applied to match the radio data), taking normalizing flux densities for the lobe regions, excluding 
the NHS and SHS, from unpublished archival low-frequency VLA (1.4 GHz, observation AL0146) 
and GMRT (610 and 244 MHz, observation number 2070) data. 3C 33 is unusual in that the flux 
density even at low frequencies is dominated by the bright hot spots, not the lobes. This calcu- 
lation predicts flux dens ities of 0.7 nJy for e ach lobe at equipartition, a factor 1.9 and 2.6 below 
the observed values. As ICroston et al.l (|2005|) find that this factor is typically ~ 2 for the sources 
they study, 3C33's lobe X-ray emission gives results consistent with those seen in other FR lis. 
The X-ray morphology does not precisely match the radio morphology, although the low X-ray 
surface brightness and heavy smoothing of the data prevent a definitive statement. Thus it is at 
least plausible that this extended X-ray emission is due to inverse-Compton emission. 



4.2.2. Thermal 

The other possibility is that this X-ray emission is from an extended hot gas corona. We con- 
sider two models for the distribution of the gas, a uniform density sphere and a shell surrounding 
the radio lobes. Modeling the emission region as a uniform-density sphere of radius 1' (68.4 kpc) 
with the best-fit temperature of the thermal model, we find a gas density of 7xl0~ 4 cm -3 and a 
total mass of 2x 10 10 M & . A more realistic density profile (i.e. a beta-model) would only change 
these numbers by a factor of order unity. This is a relatively small amount of gas, probably less than 
the mass of stars in the host galaxy. The optical luminosity of the hos t galaxy (M R =-20.7) implie s 



a stellar mass of ~9xl0 10 M assuming a mass to light ratio of 6 dBinney & Tremaine II1987|) . 
The total thermal energy and pressure of this gas is ~1.3xl0 59 ergs and ~3.2xl0~ 12 dyn cm -2 , 
respectively. The thermal energy of the gas is not partic ularly large compared with the mech anical 



energy of many powerful radio galaxies such as 3C 33 (|Birzan et al.l l2004; K raft et al.H2006|) . Any 
atmosphere of this mass, or less, is likely to have been shock heated by the radio outflow, and so 
would no longer be bound to the relatively shallow gravitational potential of the host galaxy of 3C 
33. 

Alternatively, the X-ray emission coincident with the lobes could be from a shell of gas that 
has been swept up and compressed by the inflation of the radio source. Modeling the northern 
emission region as a spherical shell with a radius of 62 kpc and thickness of 10 kpc, with the 
best-fitting temperature of the thermal model, we find a proton density in this shell of 6.0 x 10~ 3 
cm~ 3 , and a total mass of 7 x 10 8 M & . For a thinner, 1-kpc thick shell, the mass is 2 x 10 8 (similar 
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numbers are obtained for the southern lobe). As a comparison, we estimated the amount of mass 
expected to have been swept up by the expanding radio lobe, assuming the j3 model parameters 
and temperature of Section 3.4, and assuming an expanding cone of opening angle 48 degrees and 
an outer radius of 90 kpc, which gave a mass of 6 x 10 8 M . It is therefore possible that this 
emission originates from compressed gas that has been swept up by the expanding radio source; 
this gas may also have been shock heated, since the radio lobes are overpressured (see Section 4.3) 
and therefore expanding supersonically at least in their outer parts. 

We note however that if the lobe-related emission is interpreted as entirely thermal in origin, 
then the absence of inverse-Compton emission from the radio lobes would require them to be 
magnetically domina ted, which is not the case for the general population of FR II radio sources 



(ICroston et al.l 12005). We therefore argue that the majority of the lobe-related X-ray emission 
is likely to be inverse-Compton in origin; however, we cannot rule out some contribution from 
thermal emission related to radio-source environmental impact. 



4.3. Large Scale Environment and Radio Lobe Dynamics 



3C 33 appears to have a poor environment with a steep density and pressure gradient in the 
regions surrounding the radio lobes. We compared the internal minimum pressure of the radio 
lobes with the external pressure from this environment in order to investigate the dynamics of 
radio-lobe expansion. The radio lobes were modeled assuming equipartition and the same elec- 
tron spectrum and normalizing flux densities as in Section 4.2.1, and found to have minimum 
internal pressures of 9 x 1CT 13 dyn cm~ 2 (this value is slightly higher if the electron density im- 
plied by an inverse-Compton interpretation for the lobe X-ray flux is used: see Section 4.2.1). 
Our estimate of the external gas pressure depends on whether the observed diffuse X-ray emis- 
sion on large scales is thermal or IC/CMB. If we assume the emission is thermal, the X-ray data 
beyond 20" from the nucleus are consistent with a flat beta-model profile of low surface bright- 
ness. For example, a component with [3=0.5 and r c =20 kpc, normalized to the observed Chandra 
count density at 50 arcsec, would have external pressures at the midpoints of the radio lobes of 
~2xl0~ 12 dyne cm -2 , consistent with approximat e pressure balance as seen in other sources for 
which similar analyses have been performed (e.g. lHardcastle & Worrall 112000; lHardcastle et al. 

"^ 20041 among many oth- 



2002; ICroston. Birkinshaw. Hardcastle. & Worrall 



2004; 



Belsole et al. 



ers). On the other hand, if we simply extrapolate the best fit beta-model parameters derived for the 
galactic scale emission (second and third paragraphs of section 3.3), the lobes are greatly overpres- 
sured relative to the ambient medium. At the approximate midpoint of the northern lobe, a distance 



of ~ 69 kpc, the external pressure is 3.7_o 6 x 10 dyn cm (with the uncertainty dominated by 
uncertainty in the /3-model parameters); at the (projected) midpoint of the southern lobe, ~ 52 kpc 
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from the nucleus, the external pressure was found to be 7.9^} 2 2 ° x 10~ 14 dyn cm~ 2 . Although both 
lobes are likely to be in pressure balance at their inner edges, with this model for the environment, 
which has an exceptionally steep density gradient, the northern lobe appears to be overpressured 
by at least an order of magnitude halfway along its length, and both lobes must be overpressured 
at their outer edges. 



The teardrop shape of the SHS has often been interpreted as shaped by a bow shock (|Rudnick 



1988J), and supports the hypothesis that it is still moving supersonically (relative to the external 
medium) through the extended hot corona. If this interpretation is correct, the opening half-angle 
of the shock, /i is related to the Mach number of the advance of the jet head by sin(/x)~M _1 . The 
opening angle of the SHS is ~50°, implying a Mach number of ~2.4 for the advance speed of 
the jet head relative to the ambient medium. The equipa rtition pressure of the material just behind 



the radio hot spots of the SHS is ~ 10 -11 dyn cm -2 (|Rudnicklll988|) . If the bow shock inter- 
pretation is correct, this material should be in rough pressure equilibrium with the shocked IGM. 
The Mach number of the advance of the jet implies a pressure jump of a facto r of ~7 (for 7=5/3) 



acros s the bow shock in the IGM based on the Rankine-Hugoniot relations (|Landau & Lifshitz 



1989T) . The ratio of the pressures of the material behind the hot spot to the unshocked ambient 
medium is roughly consistent with this value. Assuming a gas temperature of 1 keV (consistent 
with the spectral fit above and the poor environment if all the X-ray emission is from thermal gas), 
the velocity (advance speed) of the jet head is ~1000 km s _1 . The radio outburst therefore com- 
menced ~7xl0 7 years ago if the jet has not slowed appreciably. Conservatively estimating that 
the mechanical power deposited into the gas equals the X-ray luminosity of the nucleus, the jets 
have put ~2x 10 59 ergs into the gas, a number roughly equal to the thermal energy of the gas (if 
all the diffuse X-ray emission has a thermal origin). We emphasize that this is a very conservative 
estimate: if in fact the mechanical power of the jets/lobes is 10-100 times greater tha n the nuclear 



lumin osity as has been observed in several radio galaxies in denser environments (Bir zan et al. 



2004j), the energy deposited in the gas could be considerably larger. The powerful radio outburst 



has at least doubled the energy density of any gas within ~70 kpc of 3C 33. It (the corona) is now 
likely unbound and in the process of being blown out of the galaxy. The diffuse X-ray emission 
may be the remnant of the group gas that has escaped/is escaping as a wind. 



5. Conclusions 

We have detected X-ray emission from both the north and south hot spots of the radio galaxy 
3C 33. The X-ray emission from the tip of the SHS (region SI) is consistent with the SSC pro- 
cess. For the larger regions behind SI and for the two compact regions of the NHS (Nl and 
N2), the observed X-ray flux is an order of magnitude or more above the SSC prediction. All 
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inverse-Compton scattering models are rejected unless the hot spots are far from equipartition, 
and relativistic beaming cannot be significant as X-ray emission is detected from both hot spots 
and the nucleus is heavily absorbed. The radio through optical flux density distributions show a 
gradual spectral steepening, indicating that the observed X-ray flux from regions S2-S4 cannot be 
explained by a simple power-law extrapolation of the longer wavelength spectrum. We conclude 
that the X-ray emission is most likely synchrotron, although the measured flux densities lie above 
a simple extrapolation of the radio/IR/optical flux densities, or predictions from single or contin- 
uous injection models. The distributed nature of the X-ray emission throughout the diffuse radio 
features behind the c ompact radio hot spot, not just coin cident with the radio peaks as seen in, for 



example, Cygnus A (iWilson. Young. & Shopbell Il2000|) . is analogous to the synchrotron emission 
seen in FR I jets. 

The hot spots of 3C 33 thus provide strong evidence for the presence of synchrotron and SSC 
components within the same hot spots, and for multiple synchrotron components in FR II hot spots 
where the effects of relativistic beaming are unimportant. If this radio galaxy were observed at 
lower spatial resolution and the entire volume of the SHS used for the SSC calculation, we would 
conclude that the observed X-ray emission from the SHS far exceeded the SSC prediction for the 
feature as a whole. In fact, both emission mechanisms are important, and SSC is the dominant 
mechanism in the most compact feature. This should warn us about stretching assumptions about 
uniform filling factors for the hot spots of more distant radio galaxies too far. 

High resolution UV images of the two hot spots could be made to further evaluate the conclu- 
sions of this paper. Similarly detailed, multiwavelength (IR and optical) studies of the NHS should 
also be made to confirm the synchrotron origin. Detailed optical studies of the NHS have not been 
previously undertaken probably because of the large difference in radio flux density relative to the 
SHS. Higher spatial resolution optical (i.e. HST) observations of the NHS may be particularly 
useful in constraining the morphology of this emission. 

It is now clear that the SSC explanation for the X-ray emission from the radio hot spots of FR 
lis is not the complete picture, and that synchrotron emission must be important in at least some 
of these sources. We have demonstrated that at least in the case of 3C 33, relativistic beaming 
is unlikely to play an important role. The role of beaming in general, however, is not known, 
and can only be addressed by studying statistically complete samples. The relative importance 
of synchrotron versus SSC X-ray emission in FR II hot spots could be conclusively resolved if 
a complete sample of such radio galaxies were observed by Chandra. We believe it is crucial to 
Chandra's legacy that such a sample is observed. 
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Fig. 1— Adaptively smoothed Chandra/ACIS-S image of 3C 33 in the 0.5-2.0 keV band with 1.4 
GHz radio contours (4'/0 resolution) overlaid. The positions of the active nucleus and the north 
(NHS) and south (SHS) hot spots are labeled. 
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Fig. 2.— Unsmoothed Chandra/ACIS-S image of the SHS in the 0.5-5.0 keV band with 4.9 GHz 
radio contours (1'.'37 x 1"29 resolution) overlaid. Regions SI through S4 are also shown. Contour 
levels correspond to 2.0, 3.7, 6.9, 12.9, 23.9, 44.4, 82.6, 153, 285, and 531 mJy beam- 1 . 
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Fig. 3.— Unsmoothed Chandra/ACIS-S image of the NHS in the 0.5-5.0 keV band with 4.9 GHz 
radio contours (1//37 x 1'.'29 resolution) overlaid. Regions Nl and N2 are shown, region N3 lies to 
the southwest. Contour levels are the same as Figure [2l 




Fig. 4.— WFPC2 (F702W filter) image of south hot spot with 4.9 GHz radio contours (0743 x 0"38 
resolution) overlaid. 
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Fig. 5.— SpitzerHRAC images of the SHS at 4.5 (left) and 8.0 (right) /xm with 4.9 GHz radio 
contours (1'.'37 x 1''29 resolution) overlaid. The radio contours are logarithmically spaced between 
3.0 and 400.0 mJy beam -1 . 
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Fig. 7. — Radio through X-ray flux density distributions of regions SI (left) and S2 (right). The 
continuous curve denotes the synchrotron spectrum, and the dashed and dotted lines the estimates 
of the SSC and IC/CMB X-ray fluxes, respectively. The opti cal and IR data for S 1 include flux 
densi ties taken from Meisenheimer et al. (1989) (as used by lHardcastle. Birkinshaw. & Worrall 
1998) as well as our new radio, Spitzer and HST values. 







5 



log(iVHz) 







5 



bg(i//Hz) 



Fig. 8. — Radio through X-ray spectral energy distribution of regions S3 (left) and S4 (right). The 
continuous curve denotes the synchrotron spectrum, and the dashed and dotted lines the estimates 
of the SSC and IC/CMB X-ray fluxes, respectively. 
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Fig. 9. — Radio through X-ray spectral energy distribution of regions Nl (left) and N2 (right). The 
continuous curve denotes the synchrotron spectrum, and the dashed and dotted lines the estimates 
of the SSC and IC/CMB X-ray flux, respectively. 
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Fig. 10. — Radio through X-ray spectral energy distribution of region N3. The continuous curve 
denotes the synchrotron spectrum, and the dashed and dotted lines the estimates of the SSC and 
IC/CMB X-ray flux, respectively. 
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Fig. 11. — Unsmoothed Chandra/ ACIS-S X-ray image of the 3C 33 nuclear region in the 0.3-0.7 
keV (left) and 1.5-2.0 keV (right) bands at full resolution (1 pixel = 0.492"). 



Fig. 12. — Adaptively smoothed Chandra/ ACIS-S X-ray image in the 0.5-2.0 keV band with all 
point sources other than the X-ray hot spots removed. 1 .5 GHz radio contours are overlaid. The two 
regions of diffuse X-ray emission, DN1 and DS1, associated with the low surface brightness radio 
emission between the hot spots and nucleus have been labeled. We attribute this X-ray emission to 
inverse-Compton scattering of CMB photons off the relativistic electrons. 
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Fig. 13. — Chandra radial surface brightness profile in the 0.5-2.0 keV band. The two dotted curves 
are the point source and beta-model profiles convolved with the Chandra PSF. The dot-dashed line 
is the sum of the two dotted curves. 
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Fig. 14. — Best-fit X-ray spectrum and residuals of the nucleus of 3C 33. The model (Model 1) 
shown in the left figure consists of a single heavily absorbed power law, a reflection component, a 
Gaussian line at ~6.4 keV, and a thermal (APEC) component. The model shown in the right figure 
(Model 2) consists of a single heavily absorbed power law, a second, less absorbed, power-law 
component, a Gaussian line at ~6.4 keV, and a thermal (APEC) component. The best-fit values of 
the parameters and uncertainties are contained in Tabled 
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Observation ID 


Hotspot 


Frequency 


Configuration 


Date 


Map resolution 






(GHz) 






(arcsec) 


AR130 


S 


4.9 


A 


1986 Jun 05 


0.43 x 0.38 


AR148 


N 


4.9 


B 


1986 Jul 13 


1.37 x 1.29 


AR148 


N 


15.0 


B 


1986 Jul 13 


0.45 x 0.40 


AR148 


S 


4.9 


B 


1986 Jul 13 


1.37 x 1.29 


AR148 


S 


15.0 


B 


1986 Jul 13 


0.45 x 0.40 


AR223 


s 


15.0 


B 


1990 Sep 14 


0.45 x 0.40 


AR208 


s 


15.0 


A 


1990 Apr 24 


0.19 x 0.12 



Table 1 : Details of VLA data sets retrieved from the public archive. The resolutions quoted are the 
FWHM major and minor axes of the restoring elliptical Gaussians. 
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Region 


RA 


DEC 


Dimensions 


Angle 


Rate 


Flux Density 


Nl 


01 


08:56.4 


+13 


22:22.9 


2.5" x 2.5" 


50° 


2.7±0.8 


0.27±0.08 


N2 


01 


08:56.3 


+13 


22:24.9 


2.5" x 2.5" 


50° 


2.0±0.7 


0.19±0.07 


N3 


01 


08:55.8 


+13 


22:13.7 


1 2.4" x 24.2" 


45° 


2.5±0.7 


0.24±0.07 


DN1 


01 


08:55.4 


+13 


21:13.0 


1.4'xl.l' 


85° 


11.8±1.0 


1.3±0.1 


SI 


01 


08:50.4 


+13 


18:26.3 


2.7" x 2.5" 


40° 


1.4±0.6 


0.14±0.06 


S2 


01 


08:50.5 


+13 


18:28.5 


3.0" x 3.0" 


40° 


3.1±0.9 


0.32±0.09 


S3 


01 


08:50.4 


+13 


18:33.0 


4.7" x 8.6" 


40° 


3.5±1.0 


0.35±0.10 


S4 


01 


08:50.7 


+13 


18:42.4 


5.9" x 5.9" 


40° 


4.5±1.2 


0.45±0.12 


DS1 


01 


08:52.8 


+13 


19:24.0 


1.4'xl.2' 


36° 


16.1±1.0 


1.8±0.1 



Table 2: Summary of regions used for flux density distributions. The dimensions are in arcseconds 
or arcminutes (for DN1 and DS1), the angle is the rotation of the box east of north, rate is the net 
(i.e. background subtracted) X-ray count rate in units of 10~ 4 cts s -1 in the 0.5-5.0 keV band, and 
the flux density is in nJy at 1 keV assuming a power law spectrum with photon index 2 (1.7 for 
DN1 and DS1) and Galactic absorption. All error bars are la counting statistics uncertainties. 
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VLA (mJy) 


Spitzer (/iJy) 




1.477 GHz 


4.9 GHz 


15.0 GHz 


231 GHz 


23.9 fim 


7.87 fim 


5.73 fim 


4.49 fim 


3.55 fim 


SI 




622. 


270. 


34. 


1554±119 


320±19 


207±21 


146±9 


99±6 


S2 




506. 


169. 




777±90 


183±21 


123 ±20 


67±6 


47±4 


S3 


1670. 


617. 


158. 






<82 


<68 


<20 


<11 


S4 


435. 


140. 








<45 


<77 


<9 


<8 


Nl 




64. 


24. 




196±137 


22±20 


20±18 


18±5 


12±4 


N2 




60. 


22. 




178±125 


26±22 


32±19 


9±5 


6±4 


N3 


1360. 


460. 








<66 


<115 


<13 


<12 



Table 3: Summary of radio and Spitzer/IR flux density measurements. The uncertainties on the 
Spitzer points are la. The upper limits on the Spitzer IRAC data are 3a (statistical uncertainties 
only). 
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Optical (/xJy) 


Galex (pJy) 


Chandra (nJy) 




6696 A 


6500. A 


6493.5 A 


4800 A 


4497.8 A 


2297A 


1524 A 


IkeV 


SI 


4.5 


3.9±2.8 


3.7 


1.3 


3.4±0.7 


<2.0 


<6.1 


0.14±0.06 



Table 4: Summa ry of optical, UV, and X-ray flux density measur e ments for SI. Optica l point s 
are taken from ISimkin I (119861) : ICrane. Stockton. & Saslow I (|1987|) : iMeisenheimer et al. I (|1989l) . 
Uncertainties on the X-ray data points are la counting statistics only. The GALEX points are not 
detections, only the 3a upper limits (see text for complete discussion). 
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Model 1 


Model 2 




Power-Law plus Reflection Model 


Two Power-Law Model 


Primary Power law 


Photon Index 


1.7±0.2 


1.7±0.2 


N H 


(5.4±0.8)xl0 23 cm" 2 


(3.5±0.4)xl0 23 cm" 2 


Flux 


2.4 xKT 11 ergs cm~ 2 s _1 


2.9 xl0~ n ergs cm" 2 s" 1 


Luminosity 


2.0 xlO 44 ergs s" 1 


2.5 xlO 44 ergs s" 1 


Fe K a line 


Centroid 


6.38±0.04 keV 


6.37±0.04 keV 


Flux 


2.4xl0~ 5 photons cm -2 s _1 


3.1 xl0~ 5 photons cm -2 s _1 


EW 


210±70 eV 


170±65 eV 


Thermal Component (APEC) 


k B T 


0.42±0.06 keV 


0.46±0.05 keV 


Z]Mg 


2.0 


1.0 


Zsi 


2.0 


1.0 




0.15±0.05 


1.0 


Flux 


5.9xl0~ 14 ergs cm -2 s _1 


5.6 xlO™ 14 ergs cm~ 2 s _1 


Luminosity 


5.1 xlO 41 ergs s _1 


4.8 xlO 41 ergs s" 1 




Reflection Component 


Second Power Law 


Photon Index 


1.7 (fixed) 


2.0 (fixed) 


Relative Norm. 


2.6 




N H 




1.2±0.4xl0 22 cm" 2 


Flux 


1.6 xl0~ 12 ergs cm~ 2 s _1 


8.3 x 10~ 13 ergs cm~ 2 s _1 


Luminosity 


1.2xl0 43 ergs s" 1 


7.0 xlO 42 ergs s" 1 


X 2 /dof 


27.3/18 


48.1/19 



Table 5: Summary of best fit spectral parameters of the nucleus of 3C 33. All fluxes are unabsorbed 
in the 0.1-10.0 keV band, all luminosities are in the 0.1-10.0 keV band. The equivalent width (EW) 
of the Fe K a line and the normalization of the reflection component are relative to the primary 
power law. All uncertainties are at 90% confidence for one parameter of interest. 
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Synchrotron 


ssc 


Region 


P 


AP 


^break 


Radius 


Length 


B eq (pG) 


SI 


2.5 


2.7 


2xl0 5 


Section 4.1 


195 


S2 


2.7 


3.2 


4xl0 5 


1.5" 


2.96" 


79 


S3 


2.8 


2.8 


2xl0 5 


2.35" 


8.6" 


53 


S4 


2.9 


2.8 


2xl0 5 


2.95" 


10.1" 


34 


Nl 


2.4 


1.3 


2xl0 5 


1.25" 


2.5" 


43 


N2 


2.4 


1.3 


2xl0 5 


1.25" 


2.5" 


42 


N3 


2.9 


1.3 


2xl0 5 


6.2" 


24.2" 


25 



Table 6: Summary of model paramters for synchrotron and SSC fits. Synchrotron fits assume an 
electron energy distribution of the form E~ p to the break energy, and a form E~ (yP+Aps> above the 
break. The index of the electron energy distribution, P, is determined from the observed spectral 
index of the radio flux density, a and is given by P — 2a + 1. The equipartition magnetic field, 
B eq is computed assuming no relativistic protons and unity filling factor. See section 4.1 for more 
detailed discussion of model fits. 



